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Abstract 

The LHC is foreseen to finally bring also the nuclear collisions to the TeV scale 
thereby providing new possibilities for physics studies, in particular related to the 
electro-weak sector of the Standard Model. We study here the Z and W ± production 
in proton-lead and lead-lead collisions at the LHC, concentrating on the prospects 
of testing the factorization and constraining the nuclear modifications of the parton 
distribution functions (PDFs). Especially, we find that the rapidity asymmetries in 
proton-nucleus collisions, arising from the differences in the PDFs between the colliding 
objects, provide a decisive advantage in comparison to the rapidity-symmetric nucleus- 
nucleus case. We comment on how such studies will help to improve our knowledge of 
the nuclear PDFs. 
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1 Introduction 



The production of W 1 * 1 and Z bosons are one of those processes which appear to be well- 
understood and under control within the toolkit of perturbative QCD (pQCD). The 
calculations are already reaching the NNLO accuracy p] with the NNLO correction 
on top of the NLO one found rather small. This observation has given confidence to 
include the recent Tevatron measurements [H [31 H] for W ± and Z production to the 
latest global analyses of the free proton parton distribution functions (PDFs) [HI El E] - 
These particles will also be copiously produced at the LHC, and are presumed to 
provide stringent tests of the Standard Model and eventually also new constraints to 
the PDFs. 

In contrast to the global analyses of the free proton PDFs [51 El E] the corresponding 
nuclear ones [SI El [101 ED] suffer from a much more limited amount and type of exper- 
imental data. Consequently, the tests of the factorization — a working hypothesis 
typically assumed in nuclear collisions — although in general fulfilled, are weaker. The 
present status is that one of the sets of nuclear PDFs, EPS09 [8], is able to satisfactorily 
combine DIS data on nuclear targets, Drell-Yan dilepton production in proton-nucleus 
collisions and inclusive pion production in deuteron-gold (dAu) collisions in central 
rapidities at RHIC with no large tension observed among the different data sets. In 
addition, the same set provides a good description of neutrino DIS data from nuclear 
targets |12j*l offering further strength to the hypothesis of the factorization. All anal- 
yses include DIS with nuclei and Drell-Yan production in proton-nucleus collisions, as 
already done in the pioneering global fits [T51 [IB] , while dAu data is only introduced 
in the EPS-framework, after the EPS08 set in [T7] . 

Along with the unprecedented center-of-mass energy, the LHC will require a good 
knowledge of the nPDFs in a kinematical region never explored before by DIS or other 
experiments in order to disentangle the physics signals from the backgrounds. It is 
therefore of utmost importance that checks of factorization are performed also within 
the LHC-kinematics. One of the goals of the planned proton-nucleus program at the 
LHC is precisely to provide the experimental constraints needed for this purpose. It is 
worth mentioning that the relevant kinematical region of nPDFs probed by the nuclear 
program of the RHIC-BNL was basically already covered by previous experiments, and 
therefore e.g. the discovery of the strongly interacting medium in gold-gold collisions 
there was something that could not be questioned by the insufficient knowledge of the 
nPDFs. 

In this paper, we study the prospects of observing the nuclear effects — the dif- 
ference between cross-sections in collisions involving heavy ions and those involving 
only free nucleons — in production of heavy bosons at the LHC0. We do not sepa- 
rately consider the different experiments (e.g. cuts imposed on the leptonic transverse 
momentum), or centrality dependence (see e.g. [20]), but we will compute the mini- 

1 See, however, Refs. [13] and [14] for contradicting result. 

2 Remarkably, the first Z data in nuclear collisions have already been published by the ATLAS and 
CMS collaborations [HI [19]. 
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mum bias rapidity y R spectra of the decay lepton-pairs in the invariant mass M 2 bin 
of the heavy boson peak, M 2 = M^ orW± . In reality, a detailed comparison with the 
experimental results will require an elaborate calculation with different cuts applied. 
However, the size of the expected effects should be well-estimated by the calculations 
presented here. For example, we have checked that integrating over the heavy boson 
resonance profiles or considering charged lepton production instead of lepton-pairs in 
the case of W decay would not change the main results. 

At the moment of writing this paper, the foreseen nucleon-nucleon center-of-mass 
energies yfs and luminosities C in different types of collisions at LHC — setting the 
stage for our studies reported in this paper — are the following ones: 

System y/s Luminosity C 



Proton + Proton (pp) 7TeV and 14TeV 10 34 cm" 2 s _1 
Lead + Lead (PbPb) 2.7 TeV and 5.5 TeV lO^cnrV 1 
Proton + Lead (pPb) 8.8 TeV lO^cnrV 1 



With clearly smaller nucleon-nucleon luminosities [21] and shorter envisaged running 
times (t ~ 1 month = 10 6 s), one could worry whether the integrated yields in pPb- 
and PbPb-collisions will be large enough to realize the studies presented here. In the 
case of pPb collisions the number of expected events would be 

W pPb „ / a pPb \ 208 /a pPb 

208 x £ x t x 



Ay R AM \Ay R AM ) ' 10 V P b 

Where <j pPb refers to a cross-section per nucleon — the ones we will report in the rest 
of the paper. In the narrow resonance approximation [1] , 



AA> Pb 7rT v Af pPb 

x 



Ay R 2 Ay R AM 



M=M V 



where Ty is the width of the heavy boson. As an example, for the Z-production this 
translates (see Fig. [2] ahead) to around 4000 events in the mid-rapidity y R = unit 
bin and still around 1000 events at \y R \ =3 unit bin. Moreover, knowing that larger 
luminosities could be possible [22], the heavy boson study in proton-nucleus collisions 
at the LHC looks feasible. For PbPb-collisions the corresponding numbers are similar 
as the smaller luminosity is compensated by a additional factor of 208. 

In what follows, we first describe the ingredients of our calculations, and then 
present the results for pPb and PbPb collisions for Z and W ± bosons respectively. The 
last section summarizes the main results. 
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2 Framework of the Calculations 



2.1 Z-Production 

The reaction we consider is the Drell-Yan production of lepton pairs with invariant 
mass M 2 and rapidity yn (corresponding to the invariant mass and rapidity of the 
force- carrying boson in the absence of electroweak final state radiation), in a collision 
of two nucleons Hx and H 2 at high center-of-mass energy yfs: 

d 2 a (yg. Hi + H 2 -» 1+ + t_ + X) 

dM 2 dy R ' ^ J 

The framework of the calculations is the NLO-level pQCD improved parton model. 
The relevant expressions can be found for example from [I] or obtained from [23l [24] 
with obvious change of couplings indicated below. In order to illustrate the ingredients 
and to exactly specify the computed quantity, we record here the leading order part of 
the cross-sections. For this particular process it reads 

where q^(x\,Q 2 ) and q^(x2,Q 2 ) are the PDFs of the colliding nucleons read at mo- 
mentum fractions x\ t 2 = (M/^)e VH, ~ yH , and factorization scale Q 2 . In this work we 
always choose Q 2 f = M 2 . The electroweak couplings c 2 can be expressed as 



M 2 (M 2 - Ml) 

c q - e g 2 X e ?(M2 _ M|)2 + M|rf ^ (3) 



M 4 

+ ^Tm^mWTmW z {A] + VI){AI + VS) 



X = (v / 2G F Mf)/(16 



(4) 



em/) 



where T z = 2.4952 GeV is the total width of the Z-boson, a em = 1/137.03604 the 
QED coupling constant, Gp = 1.16637 x 10~ 5 GeV -2 the Fermi constant, and Mz = 
91.1876 GeV the mass of the Z-boson. The first term in Eq. (JH]) involving only the 
quark charge squared e 2 is the pure QED part corresponding to the one virtual photon 
exchange. The third term origins correspondingly from the pure Z-exchange, and the 
second one is the Z7-interference term. The vector and axial couplings are 

V t = -l + 4sin 2 6> w , A t = -1 

V qn = l-j|sin 2 flw, A qu = l (5) 



4 

V gd = -1 + - sin 2 W , Ag d = -1, 
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for lepton £, up-type quarks q u = u, c, and down-type quarks q u = d, s, b. The sin 2 # w = 
0.23143 is the weak-mixing angle. 

Let us note, that since we are mostly interested in certain ratios of the cross-sections, 
the sensitivity to the specific choice of the numerical values of the couplings indicated 
above, is not that critical. Also, the effect of fixing the factorization scale Q 2 to the 
invariant mass of the lepton pair, does not bear large consequences to the cross-section 
ratios — especially not to the nuclear effects — as the Q 2 -dependence of the nuclear 
effects is rather mild at such large values of Q 2 . 



2.2 W ± -Production 

In the case of W ± -production, we consider the leptonic channel 



dM 2 dy R 



(6) 



assuming that the missing neutrino momentum could be reconstructed by the experi- 
ment. In reality, extracting the rapidity spectrum for is not possible without some 
theory input (2j |26] as the neutrino momentum cannot be fully reconstructed. In this 
respect the single lepton cross-section d 2 o j (dy^dp^) for Hi + H 2 — > £ ± + X should be 
more useful as a better defined quantity from the experimental point of view. 

However, at the lepton transverse momentum p l T « Mjy/2 — where the contami- 
nation of the charged lepton spectra from heavy quark decays should be negligible — it 
can be argued (see the Appendix) that the Hi + H 2 — > + X cross-section is sensitive 
to the same partonic combination as the Hi + H 2 — > £ ± + v + X process. Indeed, we 
have checked that the main results for the ^ ± production at p 1 ^ fa Mjy/2 are essen- 
tially unchanged with respect to the W ± production at M = Mw considered here. As 
demonstrated in the Appendix, even if a wider p 1 ^ -window is utilized, the main results 
are still well-estimated by the calculations presented here. 

In leading order, the cross-section for the process of Eq. (E]) reads 



d 2 a w± Tia 2 em M 4 



dM 2 dy R 36sM 2 sin 4 6 W (M 2 - M 2 ,) 2 + M^T 



2 



(7) 



x 



\Vv\ 2 ql l \x u Q 2 f )qf\x 2 ,Q 2 f ) +qf(x u Q 2 f )q^(x 2 ,Q 2 



fJ 



where VijS denote the elements of the CKM-matrix, and 



i — u,c j = d,s,b for W + 
i = d,s,b j = u,c for W~. 

We take the W mass to be M w = 80.398 GeV, the total width T w = 2.141 GeV, and 
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Figure 1: The typical shape of the nuclear modifications as a function of x, and the kine- 
matical reach for Z-production at three different center-of-mass energies and M 2 = M| for 
some values of rapidity. 

set [27] 
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2.3 Nuclear PDFs 

The usual way of constructing a set of nuclear PDFs is to take a set of free proton PDFs 
on top of which nuclear modifications are implemented. The EPS09 set jH], employed 
for the calculations in this paper and which allows for a proper error analysis, has been 
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obtained using CTEQ6.1M [25] as the reference proton distributions. Here we will use 
instead CTEQ6.6 which includes a more elaborated prescription for the treatment 
of the heavy quarksj. To be precise, the PDFs used to compute cross-sections involving 
heavy ions are 

Z N 

jA yproton,A _|_ j neutron, A 

A A 

where Z and N are the number of protons and neutrons in a nucleus with mass number 
A. The bound proton and neutron PDFs are denoted by p roto ^ A an d yneutron,A^ tlie 
latter ones obtained from those of a bound proton by a basis of the isospin symmetry, 

^neutron, A ^proton, A ^neutron, A ^proton, A (10) 

For other flavors the bound proton and bound neutron PDFs are equal. In EPS09, the 
nuclear modifications to the PDFs are defined through multiplicative, scale and flavor 
dependent factors Rf(x,Q 2 ): 

f pvoton ' A (x, Q 2 ) = R f {x, Q 2 )/ frccproton (x, Q 2 ). (11) 

As a concrete example, the down valence quark distribution dy = d — d would be 



7 N 

t A r> jfrcc proton . r> free proton / 1 o \ 

d V = -jRd V dy + -jRuyUy P • (12) 



To understand the shapes of the calculated nuclear modifications presented in 
the following sections, we plot the typical modifications of the corresponding PDFs 
Rf(x, Q 2 ) in Figure [TJ This figure also indicates the values of momentum fractions £2 
that enter to the leading order term of Z-production Eq. (JT|) (they are also the small- 
est 22-values involved in the NLO-level calculations involving integrations) at three 
different center-of-mass energies foreseen to be realized at the LHC nuclear program. 

In what follows, we will especially consider the uncertainties related to the PDFs 
and to their nuclear modifications. To this end, EPS09 and CTEQ6.6 contain error 
sets of PDFs forming an inseparable part of the published best-fits So- Here, we 
use the symmetric prescription 



AX = 2yS[ x( ^ ) - x( ^" ) ] 2 ' (13) 

to calculate the uncertainty related to the quantity X, which could be either a cross- 
section or a combination of them. If the experimental results end up being more 
accurate than the predicted uncertainty bands calculated in this way, the data will be 
evidently capable to set useful constraints to the PDFs. 

The error analysis carried out in this manner provides a propagation of the exper- 
imental uncertainties included in the corresponding global fit of PDFs. It is, however, 



3 We have checked that the change in the nuclear ratios EPS09 or the error analysis are not very 
much affected if the global fit is performed by taking CTEQ6.6 as the reference proton set. 
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known that these global fits can be biased by the assumptions in the initial conditions 
to be used in the DGLAP evolution equations. These biases can be made smaller when 
the amount of data is large. In the nuclear case the two main biasing assumptions at 
the parametrization scale Q 2 = Qq are: 

1. The flavor decomposition is, at the moment, very badly known: EPS09 assumed 
that R Uv {x,Ql) = R dv (x,Ql) and Ru{x,Ql) = R 2 (x,Ql) = Rg(x,Ql), at the 
parametrization scale. Although this feature to some extent disappears in Q 2 - 
evolution to the scales of heavy boson mases, we anticipate that relaxing this 
assumption would especially increase the uncertainties of W ± cross-sections, be- 
ing more sensitive to the flavor decomposition. 

2. At the moment, the experimentally constrained x-region is limited to i > 10 -2 . 
Below, the shape of the nuclear modifications and the size of their uncertainties 
are restricted by the assumed form of the fit function. Such a perversion makes 
the uncertainties to appear smaller than they actually are. For pPb collisions at 
^/s = 8.8 TeV this would happen right outside the midrapidity > 0, as can 
be read off from the Figure [TJ That is, data precise enough would also facilitate 
having more small- a; freedom in the fit functions thereby obtaining more realistic 
estimations for the uncertainties in nPDFs. 

In summary, although the uncertainties presented in the following are the best 
estimates available for the propagation of present experimental information to the 
unknown LHC regime, they are to be taken as lower limits of the real uncertainties. 

3 Results 

3.1 Z Production in p+Pb collisions 

The usefulness of Z-production in pPb collisions lies in the almost symmetric rapidity 
spectrum predicted in the absence of nuclear effects in PDFs. This is demonstrated 
in Figure [2} where the dashed line and the green band represent the prediction and 
its uncertainty without nuclear effects in PDFs. The small departure from a totally 
symmetric spectrum is only due to the different relative content of u and d quarks in 
the proton and the lead nuclei. The evident smallness of such isospin effect can be 
traced back to the electroweak couplings. Especially since 




(14) 



and 




0.13, when M 2 = Mf 



0.60, when M 2 < Mf , 



(16) 



(15) 
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Z-Spectrum, pPb at 4~s = 8.8TeV, M=M Z 
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Figure 2: The calculated cross-sections for Z-boson production in pPb collisions at y/s = 
8.8 TeV at the Z-pole M 2 = M|. The dashed line represents the central prediction calculated 
without applying the nuclear corrections to PDFs, and the green band is the uncertainty range 
derived from CTEQ6.6. The solid line is the prediction computed by CTEQ6.6 applying the 
nuclear effects from EPS09. The error bars quantify the uncertainty derived from EPS09 
uncertainty sets. The red dashed-dotted curve is the prediction with only QED couplings 
(first term in Eq. [3|) multiplied by 1100 with no nuclear corrections to PDFs. The lower panel 
shows the relative uncertainties with the same color codes. 

the asymmetry between forward and backward rapidities is suppressed close to the Z- 
peak with respect to photon exchange. This is demonstrated by the red dashed-dotted 
line in Figure El which represents the calculation with only QED-couplings (normalized 
by a factor of 1100 for clarity) bringing clearly out the effects associated to the different 
isospin content of the proton and the lead nucleus. 

Applying the nuclear effects to the PDFs, however, significantly modifies the almost 
symmetric distribution and the induced difference is clearly visible by eye already in 
the level of absolute spectrum. This is evident from Figure where the solid barred 
line represents the calculation with the EPS09 nuclear effects in PDFs, and the error 
bars around this line standing for its errors. 

Although the uncertainty due to free proton PDFs in the absolute spectrum in 
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Figure 3: The predicted ratio between forward and backward yields of Z bosons in pPb 
collisions at -^/s = 8.8 TeV at the Z-pole, M 2 = M\. The green solid band represents the 
prediction calculated with CTEQ6.6 without applying the nuclear effects, and the solid black 
barred line with gray shade is the prediction computed by CTEQ6.6 applying the nuclear 
effects from EPS09. The error bars quantify the uncertainties resulting from the EPS09 
uncertainty sets. 

Figure[2]is clearly non- negligible, the absence of strong isospin effects in the Z-spectrum 
in pPb collisions can be exploited to cleanly extract the nuclear corrections to PDFs. 
This can be done by simply considering the ratio between Z-production in forward and 
backward directions 

d 2 a Z,y R d 2 a Z,-y R 

dMHyJ dM 2 dy R ' ( } 

in which the uncertainties due to free proton PDFs tend to cancel to large extent, 
thanks to the symmetry property mentioned earlier. This is demonstrated by the very 
thin green band in Figure. |3J On the other hand, the asymmetry and its uncertainty 
originating from the nuclear effects is reflected as a wide barred gray band in Figure [3j 
The reason for a rather large nuclear effect is clear: at forward rapidity (yn > 0) one 
is more sensitive to the small-x shadowing (Ri < 1), whereas at backward direction 
(yji < 0) the antishadowing (Ri > 1) sets in, see Figure [TJ In the ratio of Eq. f[T7|) 
these two effects nicely add pronouncing the predicted nuclear effect. 
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W" Spectrum, pPb at = 8.8TeV, M=M W W + Spectrum, pPb at = 8.8TeV, M=M W 




Figure 4: The calculated cross-sections for W^-production in pPb collisions at y/s = 8.8 TeV 
at the W-pole, M 2 = Myy. The dashed line represents the central prediction calculated with 
CTEQ6.6 without applying the nuclear effects, and the green band is the uncertainty range 
derived from CTE6.6 PDFs. The solid line is the prediction computed by CTEQ6.6 applying 
the nuclear effects from EPS09. The error bars quantify the uncertainties resulting from the 
EPS09 uncertainty sets. The lower panels show the relative uncertainties with the same color 
codes. 

3.2 W ± Product ion in p+Pb collisions 

The extraction of the pure nuclear effects from W ± -production is not that straightfor- 
ward than isolating them from the Z-spectra. This is because the symmetry present 
in Z-spectra in the absence of nuclear effects in PDFs does not carry over to W 1 * 1 - 
production as can clearly be seen from the asymmetric green bands in Figure H] where 
the W ± -spectra in pPb-collisions are plotted. As earlier, the predictions with and 
without the nuclear effects in PDFs are plotted separately. 

In analogy to the Z-production we first consider, in Figure the ratios 

J2 W±,y R d 2 W± -y R 

dM 2 dy R ' dM 2 dy R ' 1 ' 

Due to the asymmetry present already without nuclear effects in PDFs, the green bands 
denoting the CTEQ6.6 range deviate significantly from unity and tend to be wider than 
in the case of Z-production. As the uncertainty range from EPS09 is generally larger, 
one could use these ratios — if they are measured at LHC — as an input to a nPDF 
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Figure 5: The predicted ratio between forward and backward yields of and W + (right panel) 
and and W - (left panel) bosons in pPb collisions at y/s = 8.8 TeV. The green solid band 
represents the prediction calculated with CTEQ6.6 without applying the nuclear effects, and 
the solid black barred line with gray shade is the prediction computed by CTEQ6.6 applying 
the nuclear effects from EPS09. The error bars quantify the uncertainties resulting from the 
EPS09 uncertainty sets. 
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Figure 6: The predicted ratio between forward-backward difference of W + and and W~ 
bosons (left) and the W charge asymmetry (right) in pPb collisions at y/s = 8.8 TeV. The 
green solid line represents the prediction calculated with CTEQ6.6 without applying the 
nuclear effects. The solid black barred line is the prediction computed by CTEQ6.6 applying 
the nuclear effects from EPS09. The error bars quantify the uncertainties resulting from the 
EPS09 uncertainty sets (not visible for the W charge asymmetry). 
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fit. There is, however, another combination of W 1 * 1 cross-sections that is much more 
inert to the underlying set of free proton PDFs and its uncertainties, namely 



2„W+ 



d 2 a 



UR 



2^W+, 



d 2 a 



-<JR 



dM 2 dy R dM 2 dy R 



I 



It - 



-Dr 



dM 2 dy R dM 2 dy R 



(19) 



This quantity is plotted in Figure El revealing how already at midrapidity, where the 
cross-sections should be well measurable, there should be a large effect emerging from 
nuclear modifications to the PDFs. In fact, utilizing the leading order formulas, one 
may easily show that in the absence of nuclear modifications in PDFs this quantity is 
exactly unity and the deviations from unity are solely due to the nuclear effects in the 
PDFs. 

We note here that the usual W charge asymmetry defined as 



R 



d 2 a W\y R d 2 a W-,y R 

dM 2 dy R ~ dM 2 dy R 



I 



d 2 a w+ ^ d 2 a w ~^ 
dM 2 dy R + dM 2 dy R 



(20) 



is much more sensitive to the free proton PDFs than to their nuclear modifications. This 
insensitivity is illustrated in Figure |6] for pPb collisions where the nuclear modification 
to the curve is not visible. Such insensitivity to the nuclear modifications can be partly 
related to the assumptions made in EPS09 (see Sec. l2.3p about the flavor-blindess of the 
nuclear effects at the parametrization scale. At the moment it is unfortunately difficult 
to estimate how much larger the nuclear uncertainties would grow if this assumption 
would have been released. However, in the case the nuclear uncertainty proves truly 
negligible, the quantity Aw in pPb collisions could turn out surprisingly useful for the 
studies of the free proton PDFs. 



3.3 Z and W production in Pb+Pb collisions 

Before the realization of any pPb-runs, the LHC will be operated in pp- and PbPb- 
modes. Therefore, it is relevant to check if already heavy boson production in PbPb- 
collisions would provide useful information about the nuclear effects in PDFs. In 
spite of the dense QCD-matter being created in PbPb-collisions, the leptons from the 
decays of heavy bosons should penetrate practically unaffected through this medium 
(for a short summary, see [29]) justifying the interpretation based on the pQCD parton 
model. 

First, in Figured we plot the predicted absolute rapidity spectrum for Z and W ± 
bosons in Pb+Pb collisions. The total uncertainty (EPS09 and CTEQ6.6 errors added 
in quadrature) there is typically less than 10%. This uncertainty is of the same order 
as the quoted uncertainty when the cross-section is extracted from the absolute yields 
utilizing a relation [30] 

rf 2 iV PbPb = (iVcoii) ^V PbPb 

dM 2 dy R a^ lastic dM 2 dy R [ ' 
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Figure 7: The predicted Z and W rapidity spectra in PbPb collisions at the heavy boson 
pole, M 2 = M| w± with y/s = 2.7 TeV and \fs = 5.5 TeV. The green solid band represents 
the prediction calculated with CTEQ6.6 without applying the nuclear effects, and the solid 
black barred line with gray shade is the prediction computed by CTEQ6.6 applying the 
nuclear effects from EPS09. The error bars quantify the uncertainties resulting from the 
EPS09 uncertainty sets. 



where aj^ lastic is the total inelastic nucleon-nucleon cross-section and (N co u) is the 
average number of binary nucleon-nucleon collisions extracted from a Glauber model 
Monte Carlo simulation. This indicates that the theory calculations appear already 
accurate enough to serve as an independent check of the Glauber model. 

Part of the uncertainties in Figure [7] are, however, rather related to the overall 
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Figure 8: The predicted Z and W rapidity spectra in PbPb collisions normalized to the 
integrated cross-section at the heavy boson pole, M 2 = M| w± (px = for charged 

lepton production) with yfs = 2.7 TeV and y/s = 5.5 TeV. The green solid band represents 
the prediction calculated with CTEQ6.6 without applying the nuclear effects, and the solid 
black barred line with gray shade is the prediction computed by CTEQ6.6 applying the 
nuclear effects from EPS09. The error bars quantify the uncertainties resulting from the 
EPS09 uncertainty sets. The symbol cr tot refers to the cross-section integrated over the 
rapidity. 



normalization and not so much to the shape of the predicted spectra. In this respect 
good observables to consider are the rapidity spectra of Z and W ± bosons in Pb+Pb 
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collisions normalized to the integrated cross-section, 



dMHy R ~ a t p f b (yii^) dMHy R 



(22) 



where cr P 3 bPb refers to the cross-section integrated over the rapidity. These are shown 
in Figure [SJ These quantities are good in the sense that they are free from any over- 
all normalization uncertainties i.e. can be formed equally from the absolute yields 
d 2 N PhPh /dM 2 dy R even if the nucleon-nucleon luminosity is not known accurately. As 
can be seen by comparing Figures [7J and [SJ especially the free proton uncertainties get 
suppressed in the latter one. As the nuclear modifications to the PDFs affect more the 
shape and not much the overall magnitude of the spectra, the nuclear uncertainties 
remain essentially unchanged. 

Apart from serving as tests of factorization and luminosity determination in PbPb 
collisions, the quantities plotted in Figures [7J and |8] are not well-suited to offer fur- 
ther constraints for the nuclear modifications of the PDFs: Even if the experimental 
errors were smaller than the uncertainty in Figures [7J and [HI the fact that the uncer- 
tainties stemming from the free proton PDFs and those originating from their nuclear 
modifications are of the same order makes it difficult to disentangle between the two. 

In order to further reduce the uncertainties due to the free proton PDFs and enhance 
the visibility of the nuclear modifications, reference cross-sections from pp-collisions 
would be needed. Avoiding the problems from the absolute the normalization uncer- 
tainty we look at ratios 

dM 2 dy R 1 dM 2 dy R ' 1 ' 

The problem here is, however, that the center-of-mass energies */s pp (here we use 
v /s^ p " = 7TeV) and A/spbPb in pp- and PbPb-collisions will be non-equal at LHC, and 
the uncertainties from the free proton PDFs will not therefore perfectly cancel in these 
ratios. Rather than plotting directly Eq. ( 123]) . we normalize it by the factor 

C ^ a Norm(v / ^Pp") / ^NorrrXv^PbPb) 

dM 2 dy R 1 dM 2 dy R ' 1 ' 

calculated with the central CTEQ6.6 set. This enhances the visibility of the nuclear 
effects eliminating the disbalance between the center-of-mass energies in pp and PbPb 
collisions. In other words, we calculate 

_ Gp0- N o P m(V S PbPb) ,d 2 <TNorm(\/ S PbPb) 

RpbPb = dMWn ' dMW R ' (25) 

with the relative uncertainties from Eq. (l23j) . The results obtained in this way are 
plotted in Figure El 

Interestingly, in the case of Z bosons the nuclear modifications cause the rapidity 
slope to change the sign with both center-of-mass energies. For W ± s the effects are not 
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Z Production, M=M 2 



Z Production, M=M 2 




Figure 9: The predicted ratio between the normalized Z and cross-sections in PbPb 
and pp collisions at the heavy boson pole, M 2 = M^ w± (px = Myy/2 for charged lepton 
production) with yfs = 2.7 TeV and y/s = 5.5 TeV. The green solid band represents the 
prediction calculated with CTEQ6.6 without applying the nuclear effects, and the solid black 
barred line with gray shade is the prediction computed by CTEQ6.6 applying the nuclear 
effects from EPS09. The error bars quantify the uncertainties resulting from the EPS09 
uncertainty sets. The results are obtained by first calculating Eq. ()23p and normalizing it 
then by Eq. (|24p . The symbol a tot refers to the cross-section integrated over the rapidity. 



that dramatic — the deviations from unity are mainly caused by the isospin effects. 
Although the free proton uncertainties are, on average, smaller than those from the 
EPS09, they are in various places still too large to be neglected in a parton fit. Thus, we 
are forced to conclude that — in contrast what was observed about the p+Pb collisions 
— the clean extraction of the PDF nuclear effects from Pb+Pb collisions seems difficult. 



1(3 



However, one should keep in mind that once the experimental data from the pp-runs 
are included to the free proton PDF fits, one should observe a decrease in the free 
proton uncertainties thereby increasing the usefulness of i?pbPbS presented in Figure O 



4 Conclusion 

We have reported a study of production of the heavy vector bosons (precisely to their 
leptonic decay channels) in pPb and PbPb collisions at LHC 

From the pPb-study we find that these collisions at LHC would serve as a powerful 
tool for testing the factorization and improving the present knowledge of the nuclear 
modifications to the free nucleon PDFs. Especially, there would be no need to use 
baseline pp data, but the nuclear modifications could be extracted from the pPb data 
alone with very little sensitivity to the underlying free proton PDFs utilizing the cross- 
section ratios 

u ° P Pb /" ° P Pb 



and 



dM 2 dy R 'dM 2 dy R 

d 2 a w ~>v R d 2 a w -<-y* 



dM 2 dy R dM 2 dy R 



I 



dM 2 dy R dM 2 dy R 



The former quantity can currently be estimated with an approximate uncertainty of 
10% outside the midrapidity (y R > 1), with uncertainty diminishing from that towards 
y R = 0. The latter quantity above carries around 20% uncertainty in the whole rapidity 
interval. The experimental data with errors smaller than these, would certainly offer 
a possibility to test the factorization in nuclear collisions and to constrain the nuclear 
modifications of the PDFs. It is worth stressing that these quantities do not suffer 
from any uncertainties associated to the use of the Glauber model for normalization of 
the nuclear cross-sections. 

In contrast, distilling information about the nuclear effects in PDFs from symmetric 
PbPb collisions suffers from a sensitivity to the the underlying free proton PDFs and 
their uncertainties. Moreover, the predicted uncertainties from the nuclear corrections 
to the PDFs are not very sizable. Therefore, the data from PbPb collisions will rather 
serve as a test of the factorization and check of the Glauber model in PbPb collisions 
than constraint for the nPDFs. 
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Appendix 



A relation between £ and W ± spectra 

In the leading order, the cross-section for charged lepton production can be written as 



E 



d 3 p 



sin 4 #w 



I>«i 



Pt 

dx 2 x 2 i=e 



-Vr 



X\X 2 

Qs 2 



(26) 



x- 



with 



and 



M 2 V 2 



U,C J 

d,s,b j 



\fq { t\xi)qf\x 2 ) + u 2 qf{ Xl )qf\x 2 ) 



d, s, b for £ + 
u, c for £~, 



-y/sp T xie~ 



Vr 



U 



-\fsp T x 2 e VR 



x 2 p T e 



Vr 



p T e 



(27) 
(28) 

x 2 \/s — pre VR y/s — pTe yR 

In the narrow width approximation the cross-section is dominated by the peak of the 
VT-propagator, that is, when the kinematical configuration is such that s — M^ = 0. 
At pt = Mw 1 2 one easily finds that this happens when 



M 



w 



VR 



Mi 



x 2 



w 



-VR 



u 



Ml 



(29) 



so that 



E- 



d 3 p 



oc 



qf \xi)qf\x 2 ) + qf {x 1 )qf ) {x 2 ) 



(30) 



This is the same partonic combination at the same values of x\ and x 2 that enters also 
to the leading order W production discusssed in Section 12.21 Therefore, the ratios 
of charged lepton production at y^ and px = Mw/2, are very close to those of W ± 
production at y^ and M = Mjy. Even in a wider leptonic p^- window the main results 
remain qualitatively unchanged, although the effects would naturally slightly depend 
on the actual value of the cuts. As an example, Figure [10] shows how much difference 
would be induced in a e+ ' yR /a e+ ~ yH if px > 25 GeV is used instead of px = M w /2. For 
these pT-integrated results, the MCFM code [31] has been used. 
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W* Asymmetry, pPb at s 7 = 8.8TeV 
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Figure 10: As the right panel of Figure [5l but the red lines indicating o* + > yR /a* > VR com- 
puted at NLO with the central EPS09 set and p T > 25 GeV. 
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